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I. INTRODUCTION 



The 1994 work of Seiberg and Witten [14,15] revealed the existence of a deep cor- 
respondence between supersymmetric gauge theories and integrable models [5,7,12,13]. 
However, the specific list of correspondences is still far from complete at the present time 
(c.f. [1-4] and references therein). In particular, it is still not known how to construct the 
integrable model which corresponds to the M = 2 SU{N) gauge theory with matter in the 
symmetric representation, although the spectral curve has been identified by Landsteiner 
and Lopez using M Theory [11,6]. 

The purpose of this paper is solve this problem. In [10], we had solved a similar, but 
simpler problem, which is to construct the integrable model corresponding to the SU (N) 
gauge theory with matter in the antisymmetric representation. The main new difficulty 
is the asymmetry between the orders of the zero and the pole of the eigenvalues of the 
monodromy operator at the two compactification points of the spectral curve. The desired 
integrable model turns out to be still a spin chain PmQm but whose main feature is a 
new periodicity condition linking Pn+N+2, Qn+N+2 to Pn, Qn through a twisted monodromy 
operator. Such periodicity conditions have not appeared before in the literature, and we 
take the opportunity to discuss them in some detail. An earlier proposal of how they can 
be used to construct other integrable models is in [8] , but not pursued further there. 

The mathematical problem can be formulated very simply. It is to find an integrable 
Hamiltonian system with spectral parameter x, spectral curve 

R{x, y)=y^ + fN{x)y^ + fN{-x)x^y + = (1.1) 

and symplectic form uj = Yl'i=-]^'^ ^^i^i) ^ ^('^«)" Here /iv(a^) = X^i^o'"*^* ^ generic 
polynomial of degree N, and the parameters Ui can be viewed as the moduli of the spectral 
curve. A system with the desired properties can be obtained as follows. Let QmPn be 
complex 3-dimensional (column) vectors satisfying q^Pn = and the refiexivity condition 
Pn — hp-n-i, Qn — where h is the 3x3 matrix whose only non-zero entries are 

^31 = ^22 = hi3 = 1. Let a, 6, c be 3 X 3 matrices satisfying 

0^ = 1, ab = ba, b'^ = ac+ca, bc = cb, = 0. (1-2) 



Consider the dynamical system 
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where /i^ it) is an arbitrary scalar function, and we have set m = -\- 1 for N even and 
m = — ^ + i for A?" odd. The system (1.3) appears uncoupled, but it will not be after 
imposing twisted monodromy conditions. More precisely, we have 
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Main Theorem. Let x he an external parameter, and set Ln{x) = 1 + xq^p^. Then 

(a) There are unique 3x3 matrices gn{x) = anX^ + hnX + which satisfy the periodicity 
condition 

9n+lLn+N-2 = LnQn (1-4) 

for any fixed data a^, 6^, c^., (pn, (ln)n=r~^^~^ ^'^^^ constraint q^Pn — 0. 

(b) Consider the dynamical system (1.3) with = cihjbm — bh,Cm = ch. Then the 
system is integrable in the sense that it is equivalent to the following Lax equation 

Ln = Mn+iLn - LnMn (1.5) 

where Mn{x) is the 3x3 matrix defined by 

Mn{x) = X{ Tf. 

(c) The spectral curve V = {{x,y);det(i/I — gn{x)Ln+N~3{x) ■ ■ ■L„(x)) = 0} for the Lax 
equation (1.5) is independent of n. It coincides with the Landsteiner- Lopez curve (1.1), 
and the system (1.3) is Hamiltonian with respect to the symplectic form ui on the reduced 
phase space un = un-i = 0. The Hamiltonian is H = un-2 ■ 

II. The LANDSTEINER-LOPEZ CURVE FOR 
THE SYMMETRIC REPRESENTATION 

It is convenient to list here the main geometric properties of the curve (1.1), which 
wiU henceforth be referred to as the LL (Landsteiner-Lopez) curve. It admits the foUowing 
involution ^ 

a:{x,y)^ (^-x, . (2.1) 




Points above a; = oo and x — 

Above X = oo, there are 3 distinct solutions of the LL equation, given by y ~ x^ , 
y ~ x~^'^^ and y ~ a;^. The involution a interchanges the first two, leaving the third one 
fixed. 

Above X = 0, there are also three solutions, given by y ~ x^, y ~ x^, and y ~ 1. The 
involution interchanges the last two points, while the first one is left fixed. The points 
y and y ^ x^ cross each other, but they are not branching points. 

Genus of the LL curve 

The Riemann-Hurwitz formula says that the genus g of the LL ciirve is given by 
2g — 2 = —6 + v, where u is the number of branching points (for generic moduli, the 
branching index is 2, which we assume). The branching points correspond to zeroes of 
dyR{x,y) 

Sy^ + 2yfN{x) + fN{-x)x'^ = 0. 
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To determine their number, we determine the number of poles of dyR{x,y). These occur 

ai X = oo. At X = oo, the three sohitions y ~ , y ~ x^, and y ~ x~-^~^'^ contribute 
respectively 2N, N + 2, and + 2 poles, for a total of 4N + 4 poles. Thus there are 
also 4A^ + 4 zeroes. At x = 0, there are two zeroes y ^ x'^ and y ~ x^, at each of which 
dyR{x,y) vanishes of second order. Thus the number v of branching points is given by 
v = AN + A-2-2 = AN, and consequently 

genus{T) = 2iV - 2. 

Genus of the quotient curve 

Let (7o be the genus of Fq = F/cr. Since F has two branch points over Fq, namely 
y x^ dX X — oo and y — x^ dX x — the Riemann-Hurwitz formula applies and gives 
2^ - 2 = 2(2^0 - 2) + 2, from which it follows that 

^0 = AT - 1. (2.2) 

General case 

The LL curve can be seen as a special case of a general family of curves defined by 
the equation 

R{x, y) = y^ + fN{x)y'^ + gN+2{x)y + reix) = 0, (2.3) 

where gN+2 and rQ{x) are polynomials of degree A^ + 2 and 6, respectively. This family 
has 2 A" + 11 moduli. The genus of a curve F defined by this equation can be found as 
before. It equals 

genus{T) = 2N (2.4) 

III. CONSTRUCTION OF THE SPIN CHAIN 

We give now the proof of the main theorem. Since finding the desirable integrable 
model is an essential component of our result, we construct the model gradually instead of 
proceeding from its final description. It is natural to look for a spin chain of 3-dimensional 
vectors with a period of A^ + 2 spins, in order to arrive at a spectral curve of the form 
(2.3). The difficult steps are to create an involution of the form (2.1) and to obtain the 
correct number of degrees of freedom. 

The spin chain system 

We look for a spin chain system of the form 

Ipn+l = Ln{x)lpn (3.1) 

with the operators Ln{x) given by L^i^x) = 1 +xqnPn where x is an external variable, and 
QniPn ai'e 3-dimensional complex vectors satisfying the condition ql^Pn — 0. The vectors 
qn and pn should be viewed as column vectors, so that q^Pn is a scalar, while qnP^ is a 
3x3 matrix. 
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Twisted monodromy conditions 

The key feature of the construction is the imposition of suitable twisted boundary 
conditions. Now the usual periodicity condition Ln+N-2{x) = Ln{x) can be expressed 
as TL — LT, if we define the monodromy operator to be {Tilj)n — 'il'n+N-2- For the 
Landsteiner-Lopez curve, we require a twisted periodicity condition of the form 

gn+l{x)Ln+N-2{x) = Ln{x)gn{x) (3.2) 

with the gn{xys suitable 3x3 matrices to be chosen later. This requires in turn the 
following more subtle choice of monodromy operator Tn{x) 

N-3 

Tn{x) = g„{x) JJ Lk+n{x) (3.3) 

fc=0 

where, by convention, the indices in the product of the L^+^'s are in decreasing order as 
we move from left to right. The twisted periodicity condition (3.2) is then equivalent to 

Tn+lLfi = L^Tn, (3-4) 

which implies that the eigenvalues of T„ are independent of n. We may thus define the 
spectral curve of the system L„ by 

r = {(a;,y); det{yl - T^ix)) = 0}. (3.5) 

Construction of gn{x) 

We look for gn{x) under the form gn{x) = anx'^ + b^x + Cn, where On, hn-, Cn are 3x3 
matrices. The periodicity condition gn+iLn+N-2 = -^nfi'n is equivalent to the following 
system of equations 

Cn+l — 
T T 

(ln+iqn+N-2Pn+N-2 — QnPn'^n . . 

T T 
On+l + bn+iqn+N-2Pn+N-2 = f^n + QuPn^n 

On+1 + Cn+l(ln+N-2Pn+N-2 = + QnPn^n 

We claim that this system can be solved completely in terms of the following parameters 



(3.7) 

{Pr, Qr),---, {Pr+N-3, qr+N-s), (InPn = 0, r<n<r + iV-3, 

for any choice of initial index r. To see this, define pn+N-2, qn+N-2 by 

Pn+N-2 =Pn^n, (ln+N-2 = X~^a~^qn, (3.8) 

with Xn a scalar yet to be determined. Then p^j^^_2qn+N-2 = X~^q^Pn = and or- 
thogonality is preserved. With c„ = cq for all n and pN-2+m QN-2+n defined already as 
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indicated, the last two equations in (3.6) can be viewed as recursion relations defining 
a^+i and bn+i- Our task is to show now that A„ can be chosen so as to satisfy the second 
equation in (3.6), which we rewrite as 

Now the recursive equation for a„+i implies that 

an+ldn^qn + hn+l\~^ q^p^ana'^ = gn + quPn^nttn^ qn- 

The second term on the left hand side vanishes since p^qn = 0. Furthermore, the term 
p^hna,~^ qn on the right hand side is a scalar, so that the preceding equation implies that q^ 
is an eigenvector for the operator a„_|_ia~^. Thus the second equation in (3.6) is satisfied 
by choosing A„ to be the corresponding eigenvalue 

an+ia~^qn = \nqn-, Xn = + PuhnO-n^ q-n, (3.9) 

completing the recursive construction. 

Note that the spectral curve corresponding to generic chain constructed above has the 
form (2.3). The dimension of the phase space equals D = 27^ 6{N -2) - (N -2) - (N - 
2) — 8 = +11, which is equal to the dimension of the Jacobian bundle over the family 
of curves defined by (2.3). 

Involution on the spectral curve 

We turn to the task of choosing the twisted monodromy so that the spectral curve 

admits the desired involution (2.1). Recall that the matrix h is given by hij = 0, except 
for his = ^22 = hsi = 1. In particular, = 1. Let us impose the following constraints on 
the spin chain and the twisted monodromy conditions 

Pn = hp_n-i, qn = hq_n-i, gn{-x)hg_n-N+2{x)h = x'^ . (3.10) 

The first two constraints imply 

Tn{-x) = gn{-x)h I W L_l_^_-^{x) J h = gn{-x)hT_l^_^^^{x)g_n-N+2{x)h 

\k=Q J 

Therefore the last constraint implies that the spectral curve F admits the involution 
(x,j/) — > ^— a;, Here we made use of the fact that Ln{—x) = Ln{x)~^, which fol- 

lows at once from the orthogonality condition q^Pn = 0. For generic choice of initial index 
m in (3.10) the second constraint is non-local in term of the corresponding parameters 
(3.7). It becomes local for a special choice of m. Let us assume for simplicity that N is 
even. Then we choose m — —N/2 + 1. The constraint (3.10) for n = m is the following 
equation for g^ 

gm{-x)hgm{x)h = x^ 
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Then the last matrix equation in (3.10) is equivalent to the system of equations (1.2) for 
the matrices a = a^/i, h = and c = c^/i. The last equation in (1.2) implies that c is 
a traceless rank one matrix. Hence, it can be written in the form 

c = a0^ , 0^a = 0, 

where a, (3 are orthogonal three-dimensional vectors. The third equation can be solved for 
h in the form 

h — n {aaj3^ + a0^ a) , fi^{/3'^aa) = 1 

All the equations are satisfied for any a, /3, and any choice of a such that = 1. We 
obtain in this way the crucial fact that the dimension of the admissible set of initial data 
gm ^ {0"m: bm, Cm) is equal to 8 4 + 4. The first term is the dimension of matrices a 
and the second term is the dimension of orthogonal vectors a, (3 modulo transformation 
a — > KCK, (3 — > K~^(3. 

Degrees of freedom of the system 

The system of vectors qn,Pn with the orthogonality constraint has 5(A^ — 2) degrees 
of freedom. The symmetry condition (3.10) reduces it to |(Ar — 2) (for say, N even). Now 
the system has the following gauge invariances 

• {QniPn) — {l-i'nCln-i l^'n^Pn)- This removes i(Ar — 2) degrees of freedom. 

• a global invariance {qn,Pn) — ^ iy^~^(lmPn) under 3x3 matrices W satisfying Wh = 
hW. Such matrices W are of the form 

(Wu Wi2 Wi3\ 
W21 W22 W21 
Wis Wl2 Wii J 

Their space is 5-dimensional. However, one degree of freedom has already been ac- 
counted for since diagonal matrices are of the form of the preceding gauge invariance. 

Thus the total number of gauge invariances is ^{N — 2) +4, and the number of degrees 
of freedom for the variables Pn is 2N — 8. A similar counting also produces the same 
number 2A'" — 8 of degrees of freedom for the system when N is odd. Now, as we saw in the 
previous section, the system am, bm, Cm has 8 degrees of freedom. Altogether, the number 
of degrees of freedom of our dynamical system is then 

ifdegrees{ao, bo, cq; {qo,Po), ■ ■ ■ {qN-3,PN-3)} = 2iV, (3.11) 

which is the same as the dimension of the geometric phase space constructed out of the 
curve r/o" and its Jacobian. 

The dynamical equations of motion for qn,Pn 

The equations of motion are determined by the matrix M„ completing into a Lax 
pair with equations of motion L„ = M„+iL„ — L^Mn- In this case, they are given by the 
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matrices M„(a;) in (1.6). We claim that the matrices M„ satisfy the following periodicity- 
condition 

Mn+N-2{x) = a-^Mr^{x)an. (3.12) 
In fact, the periodicity conditions for qn and pn imply that 

-1 ^ 



Mn+N-2{x) = X 



Using the fact that qn-i is an eigenvector of a^a^l^? ^^e first term on the right hand side 

T 

can be easily recognized as "^j"^^" an- Similarly, the second term can also be rewritten 



as a^^ '^T^"' ^ an, using the fact that is an eigenvector (on the left) of the matrix andn+i 

pZananXi = KPn- (3.13) 
To prove this identity, we use first the recursive relation defining an+i and obtain 

p^ttn+i + -^plbn+ia'^qnP^an = plan- 

This implies already that is an eigenvector 



Tl —1 

^anttn+l = (1 7 ) Pn 



Pn—ii—n-tL \- \ 

It remains only to simplify the expression for the eigenvalue. This is done using the 
recurrence relation defining bn+i 

K-Pnbn+ian^qn =1 + PniK - bn+l)an^qn 

= + PliT-Cn+ian^qnPnf^n - qnPnCn)an^qn = 1- 

The proof of the relation (3.13) and hence of the periodicity relations for Mn is complete. 

Equations of motion for am, bm, Cm 

Let be the solution of the equations '^n = -^n^n, dt^n = Mn^n, which is the 
eigenvector for the monodromy matrix 

y*n = gn^n+N-2- 

Taking the derivative of the last equation we obtain 

gn = MnQn " gnMn+N-2, (3.14) 
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which is equivalent to the equations 



. _ I gn-lPn 'lnPl-l \ ^ _ -l f Qn-lPn QnPl-1 
.PiQn-1 pi-lQuJ \PiQn-l pi-lQn 

Qn-lPn QnPl-A _^ ( q^-ipl QnPl-A (3.15) 



PnQn-1 Pi-lQnJ \PnQn-l Pi-lQr, 

Cn = 0. 

The key consistency condition which has to be verified is that for m = — ^ + 1, this 
dynamical system restrict to the variety of matrices a — amh, b = bmh, c — Cmh defined by 
the equation (1.2). Among these, the difficult equation to check is = ac + ca, and we 
turn to this next. Here we have assumed to be specific that N is even. The case of N odd 
is similar. 

Relations at m = — ^ + 1 



We claim that the periodicity of the system together with their involution relations 
imply the following relations at m = — ^ + 1 

'^m Qm—l — hqrn Pm'^rn ~ Pm — 1^ (o -\ c\ 

— 1 \ h \ T Tl ^^o.lOj 

'^m — '^miT'Qm—l '^m—lPm—l^rn — Pm"' 

To see this, we note that the periodicity conditions with n = m and n = m — 1 give 
respectively 

1 



Q—m — X ^rn Imj P—m — Pm'^'m 



1 



_ -1 T _ T 

Q—m—1 — , O'm—l^'m—lj P—m—1 — Pm—l^m—1 

On the other hand, the involution relation with n = —m and n — —m — 1 give 

Q—m h,QjYi — i, P—m ^Pm — 1 
Q-m-1 = hqm, P-m-1 = hpm 

Eliminating q-m, P-m, Q-m-i, P-m-i between these relations, and applying the relation 
O'n^Qn = ^nO'n+iQrii Pn^^n = ^nPn^^n+i^ wc obtain the dcsired relations. 

In terms of a, b, c, the above relations imply in particular 

Pl-ia=-^pl, a-\ra-i = qm. (3.17) 

We now claim that 

Am = Am-i = 1. (3.18) 

In fact, the first and third relations in (3.16) imply at once qm-i = Am-i(om^)^9m-i) and 
hence A^-i = 1. Next we show that Am = 1. Recalling the expression (3.9) for Am-i, we 
may also write 

Am-l = 1 + Pm-l^m-ia^-iqm-l = 1 + Pm«(^m-l^)9m- 



using the facts that p^_i = — p^a and a^"'_i9m-i = ^m-ihqm-i- We use now the 
inductive relation on the b^s 

b + ca~^qm-iPm-i(^ = (bm-ih) + qm-lPm-l^- 
Substituting in the previous formula for A^-i gives 

Xm-l = 1 + p'^a{b + Ca~^qm-lPm-l(^ - <lm-lPm-lC)qm = 1 + PmCtbqm 

since p^a and aq^ are proportional to and qm-i respectively, and Pn and g„ are 

orthogonal. Since we also know that Xm-i = 1) we deduce that pj^abqm = 0. Now the 
relation (3.9) applies to itself, giving 

Am = ^+Pmbmam(im = 1 = 1 + Pm^aqm = 1 + p'^abqm 

where we have used the equations (1.2) for a,b,c. Since p^abqm is known to vanish, it 
follows that Xm = 1- 

We can now return to the equations of motion for am, bm, Cm- The equations (3.16) 
imply 

Pm-lQm = T Pm(^qm = T PmQm-1- (3.19) 

Using (3.16) and (3.19), it is now easy to recast the equations of motion (3.15) for 
Om, bm, Cm in terms of the equations of motion (1.2) for a, 6, c 

T T 

a^Qb + bQ, b = Qc + cQ, Q ^ . 

Piqn-l Pn-lQn 

Now the compatibility condition for b is b^ = ac + ca, which implies bb + bb — dc + cd. 
Substituting in the previous formulas show that this is verified. 

IV. THE SYMPLECTIC FORM 

We turn now to the third statement in the main theorem, which concerns the Hamil- 
tonian structure of our dynamical system. Since the arguments here are very close to 
the ones in our earlier work [10], except for corrections due to the twisted monodromy 
conditions, we shall be very brief. 

As in [10], our approach is based on the universal symplectic forms obtained in [8,9] 
in terms of Lax pairs. Although our main interest is the symplectic form u defined in 
Section I, there are other symplectic forms and fiows which can be treated at the same 
stroke. Thus we define the following symplectic forms a;(£) 

1 ^ 

^{^) ^ 2 5^ Resp„0(^), (4.1) 
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where 

^{£) = (< ^n+iiQ)^Lr,ix) A SlPniQ) >k +^*k {S9k9k^) A 5iPk) —■ (4.2) 

X 

The various expressions in this equation are defined as foUows. The notation < fn >k 
stands for the sum : 

k+N-3 

<fn>k= Yl (4.3) 
n=k 

The expression ip^{Q) is the dual Baker- Akhiezer function, which is the row- vector solution 
of the equation 

rn+l{Q)Ln{z) = rn{Q), rk+N-29k\Q) = V'^kiQ) , (4.4) 

normalized by the condition 

rk{Q)MQ) = 1- (4.5) 

Note that the last term in the definition of the symplectic form reflects the twisted bound- 
ary conditions. As we shall see, that makes the form independent of the choice of the 
initial index n = k. 

We show now that the symplectic form a;(o) coincides with uj. In fact, more generally, 

2N-2 c 

The expression is a meromorphic differential on the spectral curve F. Therefore, the 
sum of its residues at the punctures Pq is equal to the opposite of the sum of the other 
residues on F. For £ < 2, the differential is regular at the points situated over x = 0, 
thanks to the normalization (4.5), which insures that Sipn{Q) = 0{x). Otherwise, it has 
poles at the poles Zi of i^niQ) and at the branch points Sj, where we have seen that V'n+i(Q) 
has poles. We analyze in turn the residues at each of these two types of poles. 

First, we consider the poles Zi of ipn{Q)- By genericity, these poles are all distinct and 
of first order, and we may write 

Sx 

ReS2,n(^) = (< Vn+l^-^nV'n >k {^QkQk^) V'fe) A —{zi). (4.7) 

X 

The key observation now is that the right hand side can be rewritten in terms of the 
monodromy matrix Tn{x). In fact, the recursive relations V'n-i-i = L^'i'n and V'n-i-i-^n = V'n 
imply that 

< tPn+l^Lnlpn >k +'^k {SgkQk^) Ipk = 

k+N-2 /k+N-2 \ / n-1 \ 

= ^k+N-2 ( n ) ^^n n + tpk+N-2 {9k^Sgk) 1pk+N-2 = 

n=k \p=n+l / \p=k j 

= ipt+N-29k^^Tkipk = S\ny. 
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In the last equality, we have used the standard formula for the variation of the eigenvalue 
of an operator, ip^ST^ipif = 'ip'^{5y)'ipk- Altogether, we have found that 

Sx 

Resz,l)(^) =dlny{zi) A — {zi). 

The second set of poles of O(^) is the set of branching points Si of the cover. Arguing 
as in [10] p. 563, we find 



Ress,n(^) = Ress, [< il^l^-^SLndipn >k +V'fc {^9k9k ^) #fc] 



dy dx 
x^dy 



Due to the identities dL{si) = dgk{si) = 0, this can be rewritten as 
ReSsifi(£) = ReSsj 



Sydx 

{i^UN-29k^^Tkdi^k) A 



(4.8) 



Next, exploiting the antisymmetry of the wedge product, we may replace 5Tk in (4.8) by 
(5Tfe — 5y). Then using the identities 

rk+N-29k\5n -5y) = d {rk+N-29k') (y - n) 

{y - Tk)dipk = {dTk - dy)ipk, 
which result from V'fc+Ar-2^fe ^(^fc - v) = {Tk - y)'^k = 0, we obtain 

Ress.l](^) = Res^. {5 {ip*k+N-29k^) {dL - dy)'ipk) A 



Arguing as before we arrive at 



dx 



The differential form 



Ress,^(£) = Ress, {i^l+N-29k ^^V'fe) A 5y—. 

X 



dx 

{^Pl+N-2Si^k) ^Sy—. 

X 



is holomorphic at a; = for < £ < 2. Therefore 

, 2N-2 , 

^Ress, {i^l^N-29k^Hk) A^y— = - ^ Res^, {i^*k+N-29k^Hk) A^y— 

X X 
Si i=l 

Using again the fact that '0Ar+fc-2fi'fc ^ ~ y~^'^ki right hand side of the last equation 
can be recognized as 



2N-2 



^ 5\ny{zi) 



A 



X^{Zi) 
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Finally we obtain 

2N 2N-2 

2a;(^) = -^Res2,l](^) - ^ReSs,J](£) = -2 ^ 5\n.y{zi) 



1=1 Si 1=1 ^ ' 



A 



The identity (4.6) is proved. 

The Hamiltonian of the Flow 

Let M.(^c) be the reduced phase space defined by the following constraints 

•^(0) = {{(ln-,Pn, an,K, Cn);UN = ac^AT-l = ttlj/G 
M(2) = {iqn,Pn;0,n,bn,Cn);Uo = Q!o, Wl = ai}/G 

where {qn,Pn, (J,n,bn, Cn) satisfy the conditions of the previous sections, G is the group of 
all allowable gauge transformations, and ckq, cki are fixed constants. 

Lemma. Let £ be either or 2. Then the equations (1.3) restricted on M.(^t) ojre Hamil- 
tonian with respect to the symplectic form a;(^) given by (4-6). The Hamiltonians Hf^i^ are 
given by 

H{0) = UN -2, H(2) = In-UAT 

Proof. By definition, a vector field dt on a symplectic manifold is Hamiltonian, if its 
contraction iQ^u)[X) — u!{X,dt) with the symplectic form is an exact one-form dH{X). 
The function H is the Hamiltonian corresponding to the vector field dt- Thus 

id^OJ^l) = ^ ^ReSp^ Ip^+iSLnlpn >k - < Ipn+iLnSlpn >k + 
a 



Equation of motion for ip^ = {Mn + lj)'4'n implies 



2_]ReSp„ < Ip^^idLntpn >fc — = < '^n+l^Ln{Mn + lA'^n >k — = 

X X 
a a 

Resp^ < ip^j^^dLnipn >k —f- 

X 

a 

We used here the equation 

^Resp^ < 'lp*^^6LnMn'lpn >fe ^ = 

X 

a 

which is valid because the corresponding differential is holomorphic everywhere except at 
the punctures. We will drop similar terms in all consequent equations. The equation of 
motion (1.5) for implies 

< V*^(^V'n >fc = < Vn+l(^n+l^n - LnMn)5ll}n >k = 

13 



= < Vn+l-^n+l'^V'n+l - IpnMnSlpn >k - < V'n+l'^n+l'^-C'nV'n >k = 

= '^/'fc+7V-2^4+^■-2'5V'fc+^f-2 - ll^kMkSll^k- < V'n+l^n+l^iv^'^n >k 

Again the last term does not contribute to the sum of residues. 
Using the equation of motion for and the equation 

ySipk = gkS'ipk+N-2 + Sgk'ipk+N-2 - Syipk 

we obtain 

V'fc {dkOk^) S^k = ^iMkSipk - yil^k+N-2Mk+N-2gk^5il)k = 

IplMkSlpk - '4>k+N-2^k+N-2Hk+N-2 - ■0fc+Ar-2-^fe+iV-2 {Ok^Sgk) 1pk+N-2 

+ipl+N_2Mk+N-2ipk+N-2Slny. 

The last term does not contribute to a sum of the residues due to the constraints S\ny = 

0{x-^) for £ = and 5lny = 0{x) for i = 2. 
The expression for iot'^ie) reduces to 

^dtUJ^^i) = l-^ReSp^ (< Vn+l^^nV'n >k +'ipk {^QkQk^) V'fe) ^^^'^ = 



a 



dx 



X 

a 

The proof can now be completed as in [10], p. 567. 

V. ^-FUNCTION SOLUTIONS 

Since the system (1.3) is completely integrable, we can obtain exact solutions in terms 
of 6'-functions associated to the spectral curve. We give these formulas here without details, 
since their derivation is entirely similar to the one in [10] pp. 557-560, taking into account 
the twisted monodromy. 

Let ipn be the Baker-Akhiezer function, which solves the simultaneous equations 
V^n+i = Lni>ni dtipn = Mnipn- Its Components ipna, 1 < q; = 3, are given by 



'ipn,ocit, Q) = (pn,ait, Q) exp [ f udflo + tdQ'^ 

\Jqc , 

eiA{Q) + tU+ + nV + Z^) e{Zo) 



eiA{Q) + Z^) e{tU+ +nV + Zq] 



Here 0{Z) if the Riemann-theta function associated to the period matrix of the spectral 
curve; A{Q) is the Abel map; V, are the vectors of -B-periods of the meromorphic 
differentials dQ,Q, dO,'^ defined by the following requirements. The differentials dQ,Q and 
dVL^ have zero ^-periods, they are holomorphic outside the two points Pi,Pz above oo 



14 



interchanged by the involution cr, with dQo having simple poles and residues ±1, while dQ'^ 
is of the form dQ~^ = ±dx{l + 0(x~^)) at these two points. The VctiQfj) are meromorphic 
functions satisfying the normalization condition rctiQ/s) = S^p and the condition that their 
divisor of poles Zq correspond to the inital data Qn(0),Pn(0) of the dynamical system. Let 
P2 be the point above 00 fixed by the involution, and let dQi be the meromorphic form 
satisfying dfli + dQ^ = dQ~^ , where (iOJ is the image of (iOi under the involution a. 

The Laurent expansion of the last factor as Q ^ -Pi defines constants Via,Wia,wh.ich 
depend only on the curve 

V2a = / d^o, Via = hnia;->Pi ( / d^o ^llix) i = 1,3 

JQa ''Qoc 

W2a = / dill, Wia = lim^^^p. ( / dfli ^x), 1=1,3 
Let $n ^ {t) be vectors with coordinates 



Then the vector pn of the spin chain is the unique (up to multiplication; different choices 
lead to diff'erent gauge choices i'n{t) in our dynamical system (1.3)) three-dimensional 
vector that is orthogonal to ^ = 2, 3, i.e. 

and the vector qn is given by the formula 



fn n— 1 



The leading coefficients of the expansion of the Baker- Akhiezer function provide also 
the expression for the variables a„. In the normalization c = (cy) with C13 = 1 and Cy = 
for all other we find 

an = $iv+n-2$n^ 

where is the (3 x 3) matrix with columns $n • 
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